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The atmospheric Cherenkov gamma-ray telescope 
MAGIC, designed for a low-energy threshold, has detected 
very-high-energy gamma rays from a giant flare of the 
distant Quasi-Stellar Radio Source (in short: radio quasar) 
3C 279, at a distance of more than 5 bilUon Ught-years (a 
redshift of 0.536). No quasar has been observed previously 
in very-high-energy gamma radiation, and this is also the 
most distant object detected emitting gamma rays above 50 
gigaelectron volts. Since high-energy gamma rays may be 
stopped by interacting with the diffuse background light 
in the universe, the observations by MAGIC imply a low 
amount for such light, consistent with that known from 
galaxy counts. 

Ground-based gamma-ray telescopes are sensitive to the 
Cherenkov light emitted by the electromagnetic showers that 
are produced by gamma rays interacting in the atmosphere. 
These telescopes have discovered, since the first detection (in 
1989) of gamma rays in this energy range (from 100 GeV to 
several TeV), more than 20 blazars, which are thought to be 
powered by accretion of matter onto supermassive black holes 
residing in the centers of galaxies, and ejecting relativistic jets 
at small angles to the line of sight (|7]l. Most of these ob- 
jects are of the BL Lac type, with weak or no optical emission 
lines. Quasar 3C 279 shows optical emission lines that allow 
a good redshift determination. Satellite observations with the 
Energetic Gamma Ray Experiment Telescope (EGRET) aboard 
the Compton Gamma Ray Observatory (CGRO) had measured 
gamma rays from 3C 279 dU and other quasars, but only up to 
the energies of a few GeV, the limit of the detector's sensitivity. 
An upper limit for the flux of very-high-energy (VHE) gamma 
rays was quoted in (OJ. 

Using MAGIC, the world's largest single-dish gamma-ray 
telescope (4) on the Canary island of La Palma (2200 m above 
sea level, 28.4°N, 17.54°W), we detected gamma rays at ener- 
gies from 80 to >300 GeV, emanating from 3C 279 at a red- 
shift of 0.536, which corresponds to a light-travel time of 5.3 
billion years. No object has been seen before in this range of 
VHE gamma-ray energies at such a distance [the highest red- 
shift previously observed was 0.212 {S}], and no quasar has 
been previously identified in this range of gamma-ray energies. 



The detection of 3C 279 is important, because gamma rays 
at very high energies from distant sources are expected to be 
strongly attenuated in intergalactic space by the possible inter- 
action with low-energy photons (j + y — > e*' + e~). These 
photons [extragalactic background light (EBL) (6)] have been 
radiated by stars and galaxies in the course of cosmic history. 
Their collective spectrum has evolved over time and is a func- 
tion of distance. For 3C 279, the range of newly probed EBL 
wavelengths lies between 0.2 and 0.8 yum (ultraviolet/optical). 
Existing instruments that are sensitive only to higher gamma- 
ray energies have so far been unable to probe this domain; by 
contrast, MAGIC is specifically designed to reach the lowest- 
energy threshold among ground-based detectors. 
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Fig. 1. Light curves. MAGIC (top) and optical R-band data (bottom) 
obtained for 3C 279 from February to IVIarch 2006. The long-term 
baseline for the optical flux is at 3 mJy. 

In observations of 3C 279 over ten nights between late Jan- 
uary and April 2006 (total of 9.7 hours), the gamma-ray source 
was clearly detected (at >6 SDs) on the night of 23 February, 
and may also have been detected the night before (Fig. 1). As 
determined by the x^ test, the probability that the gamma-ray 
flux on all 10 nights was zero is 2.3 x 10"^, corresponding to 
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Fig. 2. Spectrum of 3C 279 measured by MAGIC. The grey area 
includes the combined statistical (Icr) and systematic errors, and 
underlines the marginal significance of detections at high energy 
The dotted line shows compatibility of the measured spectrum with 
a power law of photon index a = 4.1. The blue and red triangles 
are measurements corrected on the basis of the two models for the 
EBL density, discussed in the text. 

5.04cr in a Gaussian distribution [see (7)]. Simultaneous optical 
R-band observations, by the Tuorla Observatory Blazar Moni- 
toring Program with the 1.03 m telescope at the Tuorla Obser- 
vatory, Finland, and the 35 cm Kungliga Vetenskapsakademien 
(Royal Swedish Academy of Sciences) telescope on La Palma, 
revealed that during the MAGIC observations, the gamma-ray 
source was in a generally high optical state, a factor of 2 above 
the long-term baseline flux, but with no indication of short 
time-scale variability at visible wavelengths. The observed 
VHE spectrum (Fig. 2) can be described by a power law with a 
differential photon spectral index of a = 4.1 + 0.7stat ± 0.2syst. 
The measured integrated flux above 100 GeV on 23 February 
is (5.15 + 0.82stat ± 1.5syst) x 10"'° photons cm"^ s"'. 

The EBL influences the observed spectrum and flux, result- 
ing in an exponential decrease with energy and a cutoff in the 
gamma-ray spectrum. Several models have been proposed for 
the EBL (^. All have limited predictive power for the EBL 
density, particularly as a function of time, because many de- 
tails of star and galaxy evolution remain uncertain. We illus- 
trate the uncertainty in the EBL by using two extreme mod- 
els: a model by Primack et al. ^, close to the lowest pos- 
sible attenuation level consistent with lower EBL limit from 
galaxy counts {.9..10i: and a "fast-evolution" model by Stecker 
et al. (177b . corresponding to the highest attenuation of all the 
models. We refer to these models as "low" and "high," respec- 
tively. The measured spectra of 3C 279, corrected for absorp- 
tion according to these two models, are shown in Fig. 2. They 
represent the range for the possible intrinsic gamma-ray flux of 
the source. 

A power-law fit to the EBL-corrected points (I72l l results in 
an intrinsic photon index of o-* = 2.9 + 0.9stat ± 0.5syst (low) and 
a* = 0.5 + 1.2stat ± 0.5syst (high). The systematic error is de- 
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Fig. 3. The gamma-ray horizon. The redshift region over which 
the gamma-ray horizon can be constrained by observations has 
been extended up to z=0.536. The prediction range of EBL models 
is illustrated by \8\ (thick solid black line) and [11} (dashed-dotted 
blue line). The tuned model of tl4\ (dashed blue line) represents 
an upper EBL limit based on our 3C 279 data, obtained on the 
assumption that the intrinsic photon index is > 1.5 (red arrow). Lim- 
its obtained for other sources are shown by black arrows, most of 
which lie very close to the model [14). The narrow blue band is 
the region allowed between this model and a maximum possible 
transparency (i.e., minimum EBL level) given by \8\, which is nearly 
coincident with galaxy counts. The gray area indicates an optical 
r > 1, i.e., the flux of gamma rays is strongly suppressed. To il- 
lustrate the strength of the attenuation in this area, we also show 
energies for t = 2 and t = 5 (thin black lines), again with t^ as 
model. 

termined by shifting the absolute energy scale by the estimated 
energy error of 20% and recalculating the intrinsic spectrum. 
Further discussion of the intrinsic spectrum and the spectral en- 
ergy density can be found in d?). 

The measured spectrum of 3C 279 permits a test of the trans- 
parency of the universe to gamma rays. The distance at which 
the flux of photons of a given energy is attenuated by a factor e 
(i.e., the path corresponding to an optical depth r = 1) is called 
the gamma-ray horizon and is commonly expressed as a func- 
tion of the redshift parameter ( 7i ); we show this energy/redshift 
relation in Fig. 3. In the context of Fig. 3, we make use of a 
model based on il4i with parameters adapted to the limits given 
by (75) and fine-tuned such that for 3C 279, the intrinsic pho- 
ton index is a* - 1.5. The tuning allows for the statistical and 
systematic errors (1 SD, added linearly). Although the intrinsic 
spectrum emitted by 3C 279 is unknown, a* - 1.5 is the lowest 
value given for EGRET sources (not affected by the EBL) and 
all spectra measured by gamma-ray telescopes so far ( [76b . so 
we assume this to be the hardest acceptable spectrum. The re- 
gion allowed between the maximum EBL determined using the 
above procedure and that from galaxy counts ([Sb is very small. 

The results support, at higher redshift, the conclusion drawn 
from earlier measurements (i75i) that the observations of the 
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Hubble Space Telescope and Spitzer correctly estimate most 
of the light sources in the universe. The derived limits are con- 
sistent with the EBL evolution corresponding to a maximum 
star- formation rate at z > 1, as suggested by {8j and similar 
models. 

The emission mechanism responsible for the observed VHE 
radiation remains uncertain. Leptonic emission models (as- 
suming relativistic electrons in the jet as source of the gamma 
rays), generally successful in describing blazar data [e.g., (03], 
can, with some assumptions, also accommodate the MAGIC 
spectrum. Hadronic models [involving relativistic protons, e.g. 
dTSjl] provide a possible alternative. However, a genuine test 
of the models can be only obtained with simultaneous obser- 
vations at different wavelengths, which are not available for 
the observations described here. Future tests of these models 
should use observations from sources at all wavelengths from 
radio to VHE gamma rays. In the domain of VHE gamma 
rays, we can expect important new insights by simultaneous 
observations with the Large Area Telescope (LAT), the high- 
energy gamma-ray instrument on the Gamma Ray Large Area 
Space Telescope [GLAST ([79l l1. Our observations of this dis- 
tant source in VHE gamma rays demonstrate that a large frac- 
tion of the universe is accessible to VHE astronomy. 
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Comments to the Introduction 

Ground-based gamma ray astronomy, pioneered by the Whip- 
ple Collaboration, has led to a number of important discoveries 
in recent years, showing an unexpected wealth of sources in 
the very high-energy (VHE) sky above 100 GeV. Among them, 
some twenty extragalactic sources were found, all but one are 
blazars. Blazars are now known in a broad range of luminosi- 
ties, from low-luminosity BL-Lac objects to high-luminosity 
quasars. They emit non-thermal radiation from radio waves 
to gamma rays, originating presumably in synchrotron emis- 
sion and inverse-Compton scattering. Most blazars so far de- 
tected with ground-based telescopes belong to the rare sub- 
species of BL-Lac objects, that can be interpreted as a kind 
of dying quasar running out of fuel, and thinning out at large 
cosmological distances (7 ). 3C 279 is the first blazar of the 
quasar type observed to emit gamma rays in the VHE range. It 
had been discovered as a source of lower energy gammas by 
the EGRET detector (T), with a peak flux higher than any other 
source seen by EGRET. Subsequently 3C 279 became a much- 
studied source with multi- wave length observation campaigns 
published (EHEH^HH. 

MAGIC analysis 

The observations were carried out in the "on/off mode", which 
includes "on" events taken with the telescope pointing at the 
source, and a similar amount of "off'" events taken off'-axis at a 
sky location near the source, with very similar operating con- 
ditions. The "off" events are necessary for estimating reliably 
the background of hadronic events. 

The raw data were calibrated {^ and analyzed with the stan- 
dard MAGIC analysis and reconstruction software (|7Q] i. Data 
runs with anomalous trigger rates due to bad observation condi- 
tions were rejected. The remaining "on" data correspond to 9.7 
hours, with 4.9 hours of "off" data. The standard analysis re- 
duces the image to a single cluster of pixels by removing image 
noise and tracks far from the main shower, and calculates im- 
age parameters from the cleaned image (17Tb . These are used by 
a multi-variate method to discriminate gamma rays from back- 
ground hadrons; the method is called "Random Forest" (UZl l 
and is a multi-tree classifier algorithm which calculates, from a 
combination of image parameters, a separation variable called 
hadronness. 

All image parameters were checked for consistency between 
"on" and "off" data. The excess events were obtained by sub- 
tracting suitably normalized "off" data from the "on" data, as 
shown in Fig. S 1 . The variable Alpha is one of the image pa- 
rameters. It describes the direction of the main alignment axis 
of the pixels in the image, a zero value indicating a shower 
seemingly coming (in the camera plane) from the known posi- 
tion of the source. Alpha was not included in the multivariate 
gamma-hadron separation. Cuts in hadronness and Alpha were 



optimized using data samples from observations of the Crab 
Nebula at comparable zenith angles. For the detection (and 
the lightcurve) we require hadronness<0.12 and Alpha<12°, 
resulting in an efficiency of 40% for gammas. For these cuts 
and with a minimal Size of 100 (Size is the sum of signals, in 
photo-electrons), we obtain a peak in the distribution of esti- 
mated energy at 170 GeV. 

The significance of the signal of 23 February was calculated 
according to equation 17 in (75), resulting in 6.15cr, taking the 
difference between a smoothing parabola fit to normalized "off" 
data and the "on" data. The probability that we deal with a sta- 
tistical fluctuation in ten independent observations (trial factor 
corresponding to 10 nights of observation) is 3.87x10"^, which 
corresponds to a significance of 5.11 cr. If "off" and "on" data 
are compared bin by bin without smoothing, the significance 
is 5.5cr for 23 February, and 5.1cr after correcting for trials. 
We also applied a x^ test to the results of ten nights, giving 
50.35 for 10 degrees of freedom, or a probability of 2.3x10"^, 
which in a Gaussian distribution corresponds to 5.04cr. The 
long-term detector stability has been studied by monitoring the 
Crab Nebula over a large range of zenith angles, and for sev- 
eral years. The observed integrated flux above 100 GeV has 
shown no evidence of any variations beyond what is expected 
from statistics. We can thus conclude that the observed signal 
is indeed an isolated flare from 3C 279. 



20 lb. 



(0 

§180 

O160 

140 

1 

100 
80 
60 
40 
20 



7.^^2 



H^ 



Observation time: 82 mins 
Excess events: 166 + 28 
Excess significance: 6.15 a 



% 10 20 30 40 50 60 70 80 90 

lALPHAI [deg] 

Fig. SI. Histogram of the image parameter Alpha. The image pa- 
rameter Alpha obtained for 3C 279 in the night of 23 February 2006, 
for "on" (red dots) and "off" (crosses and filled area) pointings, re- 
spectively. The dotted line is a simple parabolic fit (without linear 
term) to the "off" events between 0° and 80°, serving to smoothen 
the background for better extrapolation towards Alpha=0; the ex- 
cess events are obtained with respect to this line. 

The reconstruction of gamma ray energy uses an approxi- 
mation derived from image parameters. It is obtained by the 
the Random Forest method using Monte Carlo gamma sam- 
ples. For the reconstruction of the energy spectrum, we used a 
low-energy analysis, making also use of the timing properties 
of the air showers with a lower energy threshold (^110 GeV). 
The cuts for obtaining the spectrum were slightly loosened, re- 
sulting in a gamma efficiency of 50%. The obtained spectrum 
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was unfolded for detector effects using ( [7?b . The spectrum was 
cross-checked with two alternative methods (one of them is the 
standard analysis and a second one also uses of the timing in- 
formation). All three spectra agree well, as shown in the figure 
(main paper) as grey area, the envelope shape for 1 cr errors of 
the spectra. The predicted energy resolution is !^23%, reason- 
ably constant above 150 GeV. Systematic errors of our detector 
(from atmospheric density and transmission uncertainties, the 
photon detection efficiency, non-linearity in the signal chain, 
limits in Monte Carlo simulations, choices in analysis cuts, etc.) 
are discussed in detail in ( 175b . Uncertainties result in a possible 
systematic energy shift of up to ±20%, which translates into a 
shift of +25% in the limit set for the gamma ray horizon (en- 
ergy for which the optical depth t = 1). The uncertainty in the 
spectral slope coefficient is ±0.2. 

Attenuation of gamma rays from distant sources and EBL 
models 

Absorption of high-energy gammas by the extragalactic back- 
ground light (EBL) depends on the column density of low- 
energy photons traversed by the gamma rays. For gamma 
rays of observed energy Ey, the maximal cross section for 
pair production at redshift z occurs at wavelength A - 
1.24jum(£'y/lTeV)(l + zf. Folded with the shape of the 
cross section, MAGIC's energy range specifically probes (with 
3C 279) background wavelengths of 0.2 yum - 2 //m. This high- 
lights the importance of the ultraviolet-to-near-infrared wave- 
length range in modifying the spectra of sources at gigaparsec 
distances (76"). Knowledge of the amount of EBL along the line 
of sight to gamma ray sources is based on bolometric measure- 
ments and galaxy counts. Bolometric measurements are gener- 
ally hampered by foreground emission in our own Galaxy {T7\. 
Since the light emitted from galaxies changes with time follow- 
ing the star formation history ( [7^ and with the ageing stellar 
populations, the EBL has changed over time, too. An in-depth 
review can be found in (79). 

Various models for the EBL have been published. In (l20l l. 
several parameterized models using different assumptions for 
the evolving EBL are introduced, taking into account magnetic 
fields, gas and dust, the relative abundance of elements heavier 
than Helium, and the star formation rate. A range of parame- 
ters is proposed such that the results agree with what is known 
from deep galaxy surveys and infrared background measure- 
ments available in 2002. In a recent publication (2T), two dif- 
ferent backward-evolution EBL models were used to calculate 
the intergalactic photon density from z = to z = 6 (called 
"Base" and 'Fast" evolution). In this report we use the "Fast" 
evolution model as predicting maximally high absorption. In a 
different publication (|22l l. a semi-analytical forward-evolution 
approach was used to obtain the evolving EBL density. This 
model predicts a rather low EBL level at z = 0, just on top of the 
galaxy counts obtained from the Hubble Space Telescope dJSt . 
In fact, the model (TT) predicts the EBL level at wavelengths 
above 5 /^m slightly below the galaxy counts obtained from the 
Spitzer (|25]|2l|25ll data and the ISOCAM results ^\27^. The 
model (2T) also lies below the FIRAS measurements in the far- 
infrared, an undershoot attributed to the abundance of Ultra- 



Luminous Infrared Red Galaxies (ULIRG), only recently dis- 
covered and not taken into account in the model. However, the 
EBL domain above 3 //m does not influence our result, given 
that the highest energy measured by MAGIC in the 3C 279 
spectrum is below 500 GeV; EBL photons with A > 2.5 jum do 
not reach the energy necessary for pair production. We have, 
therefore, chosen the model {22) as the low limit for the EBL 
level in the frequency range relevant for the MAGIC 3C 279 
measurement. 

In the discussion of the gamma ray horizon in the main pa- 
per, we have used a model tuned to give for 3C 279 the intrin- 
sic photon index a* = L5. We call this ad-hoc EBL model 
maxEBL and give more details in the following. It is based 
on ( 12Q|) . with modified parameters. The "best-fit" model, de- 
fined in (28), is changed by reducing the optical and infrared 
star formation rate, by adding a "warm dust" component and by 
taking into account a high fraction of UV emission, which fits 
the data derived from AGN absorption lines, a method called 
proximity effect [e.g., {29)]. Additionally, the dust parameter 
for the optical galaxy component is set to a very low value of 
E(B-V) = 0.01, to account for the maximum possible optical 
flux. In detail, the parameters used [definition in {20}] are the 
following: 

SFRoPT : a = 3.5, /3 = -l.2,Zp = 1.2,p.(z^) = 0.08, 
SFRlig : a = 4.5,^6 = 0.0,Z;, = 1.0,p,(Z;,) = 0.09, 
and furflier/esc = 4,C2 = lQ-^^\ E(B - y)oPT = 0.01. Con- 
trary to the "best-fit" EBL, the maxEBL model is above the 
emissivity data at wavelengths between 0.16 - 1.0 jum [ (l20b 
and references therein]. 

All EBL models mentioned contain the evolution of the ra- 
diation contributing to the EBL at different redshifts, and use 
the Lambda-CDM Universe with D.^ - 0.7 and Qm - 0.3, 
confirmed by the latest results from the Wilkinson Microwave 
Anisotropy Probe (WMAP) {30). The models are graphically 
summarized in Fig. S2, along with a number of measurements. 
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Fig. S2. EBL Models. Some of the models for EBL, for z=0, and 
measurements at various wavelengths. The green line is based 
on (20j, with parameters adapted to the results from the latest 
galaxy counts. The dotted line is from J57]. The shaded verti- 
cal band indicates the range of frequencies corresponding to the 
IVIAGIC observations. 
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Comparing the effect of different EBL models on our mea- 
sured spectrum we add the following comments: 

• The high model dZTT l leads to an intrinsic spectrum of 
3C 279 with a (fitted) index of a* = 0.5 + 1.2 (see main 
paper), a rise in the VHE gamma ray spectrum very diffi- 
cult to reconcile with extrapolation from EGRET data and 
with general constraints in the spectral energy distribution. 

• The maxEBL model has been derived to be in agreement 
with an intrinsic blazar spectral index of a* = 1.5, and 
can be considered as an upper limit on the EBL density 
respecting the constraints from the measured VHE blazar 
spectra (16','31). Note that the values for the optical and 
infrared star formation rate (SFR) and for the optical dust 
parameter are close to their lower limits, and no AGN con- 
tribution has been taken into account. 

• The low EBL derived from {22) apparently gives an ac- 
ceptable result, resulting in a most probable intrinsic 
power law with a spectral index a* - 2.9 + 0.9 (see main 
paper). 

Some of the models discussed above seem less likely in the 
light of the present measurements, requiring the sources to 
exhibit a redshift-dependent hardening of their spectra. We 
should also note that the EBL density may be higher than de- 
rived here, if VHE photons are speculated to undergo interac- 
tions beyond standard physics [e.g., G52)1. 

Spectral Energy Distribution 

Blazars represent the most violent variant of the AGN phe- 
nomenon, related to the activity of supermassive (M = 10^-10^ 
solar masses) black holes harbored in the center of galaxies. 
Their emission, extending from radio wavelengths to high- 
energy gamma rays, is produced within a jet of plasma ejected 
at relativistic speeds and pointing close to the Earth. 

The spectral energy distribution (SED) of blazars is charac- 
terized by two broad bumps, the first one peaking between the 
infrared and the X-ray band, the second one in the gamma ray 
domain. The first peak traces the synchrotron emission of rela- 
tivistic electrons spiraling along the lines of the magnetic field 
in the jets; the high-energy peak is presumably due to the in- 
verse Compton scattering between these relativistic electrons 
and low-energy photons. The latter can be the synchrotron 
photons themselves [Synchrotron-Self-Compton(SSC) mecha- 
nism i33[ ] or ambient radiation [External-Compton (EC) mech- 
anism ( |54]|55l )1. Ambient radiation can enter the jet either from 
an accretion disk or from surrounding gas clouds. These mod- 
els ultimately rely on electrons (or pairs) accelerated within the 
flow and thus are called leptonic models. 

Other possibilities (hadronic models) involve gamma ray 
emission due to accelerated protons and ions. Such emission 
can result from pion production in hadronic interactions, pair 
production, or be due to synchrotron radiation [e.g., ( |36]|77||55I 
I39l l1. and includes the associated production of neutrinos and 
ultrahigh-energy cosmic rays (4U). Although the energy losses 
of protons and ions are suppressed by their large mass and the 



low matter density in the acceleration zone, they become size- 
able at ultrahigh energies where photo-production of pions can 
take place in collisions with low-energy synchrotron photons. 

Until now the majority of known extragalactic sources of 
VHE photons belongs to the class of BL-Lac objects, whose 
spectra show a weak signature of thermal radiation from out- 
side the jet and hence favor the SSC scenario. In quasars 
[mostly Flat Spectrum Radio Quasars (FSRQs)], such as 
3C 279 the situation with thermal radiation is different (the blue 
bump in the spectrum), and the EC process is likely to dominate 
if the gamma ray emission zone is close to the origin of the jet. 
The validity of the EC mechanism is not undisputed, however, 
since the gamma ray emission zone could well be far away from 
the central source of thermal radiation [e.g., ( I47ll42l i1. 

Hadronic emission models require much larger magnetic 
field strengths than leptonic models, in order to confine the 
more massive and energetic particles to the jet. They gener- 
ically produce harder VHE spectra than leptonic models. In 
low-peaked blazars at high redshifts such as FSRQs, however, 
the spectrum is steepened due to attenuation by the EBL, ren- 
dering differences in the spectra difficult to detect. 3C 279 is 
among the nearest FSRQs, and the energy threshold of MAGIC 
is low enough to permit observing this source. 

The EGRET gamma ray telescope (operating in the past 
decade and sensitive between 20 MeV and few GeV) observed 
at several occasions bright and hard (photon index around 2) 
emission from 3C 279. In particular, during an extensive cam- 
paign held in 1996 the source showed a large (a factor 10 in 
amplitude) and short (~ 2-3 days) gamma ray flare (01 with 
episodes of relatively fast variability, including an increase by 
a factor of 2.6 in ~8 hours. Naively extrapolating this behavior 
to the MAGIC band one could expect bright and short episodes 
also in this band. However, in the standard leptonic framework 
the detection of photons with energy above 100 GeV emitted 
by 3C 279 is somewhat surprising. Indeed, a careful consid- 
eration of the conditions expected in the source raises several 
problems, related to the high energy required for the electrons 
to emit photons with such high energy and the strong opacity 
to gamma rays expected inside these sources. The detection of 
strong VHE emission from 3C 279 thus is stretching leptonic 
models to or beyond their limits. 

Emission model 

The simplest EC model ("one-zone") assumes that the bulk of 
the emission is produced within a single region. In order to 
properly model the emission, good coverage of the entire SED 
is required, since only with known position and luminosity of 
both peaks is it possible to fully constrain the physical param- 
eters [e.g., (l?Jlll. In the case of the 3C 279 flare, we only have 
the simultaneous measurement in the optical R-band with the 
VHE spectrum derived by MAGIC. With this input we tried to 
reproduce the observed VHE spectrum assuming a simple one- 
zone model [see (W) for a full description]. 

The MAGIC observations contribute data to the important 
energy range above lO^^'Hz, poorly studied thus far because ly- 
ing above the range of satellite-borne instruments and below 
that of most ground-based Cherenkov telescopes. The SED in 
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Fig. S3 shows the MAGIC points corrected according to the 
two extreme EBL models discussed in the paper, together with 
measurements at lower frequencies from different epochs. In 
the same figure we show the generic one-zone model described 
here. The model uses the following assumptions and param- 
eters: The source is spherical with radius R - 5 x 10'^ cm, 
and it is in motion with a bulk Lorentz factor F = 20 at an an- 
gle 6 = F ' - 2.9 deg with respect to the line of sight. The 
magnetic field has an assumed intensity B = 0.15 G. The par- 
ticles, with a total density n = 2 x 10'* cm"^, follow an energy 
distribution represented by a broken power law extending from 
yi = 1 to 72 = 3 X 10^, with indices n\ -2 and ni = 3.7 below 
and above the break at yy, = 2.5 x 10^. We model the external 
radiation field assuming that a fraction t - 0.005 of the disk 
emission, a black body with luminosity 6x10'*^ erg/s, is di- 
luted into the Broad Line Region, a sphere with radius 4x10'^ 
cm. 
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Fig. S3. Spectral energy distribution for 3C 279. Observations 
af different source intensities, over the years from 1991 to 2003, 
from (45\ and (4S\, with IVIAGIC points (2006) at the high-energy 
end (right). The IVIAGIC points are corrected according to the two 
extreme EBL models of j22) and (27) (blue and red triangles, re- 
spectively); the former clearly gives a more probable result with a 
spectral index of a' = 2.9 + 0.9. An emission model and its indi- 
vidual components are also shown: the solid line represents the 
total model emission, fitted only to the (blue) IVIAGIC points and the 
(black) triangle from the KVA telescope. The individual components 
of the emission model are synchrotron radiation (dotted line), disk 
emission (long-dashed), synchrotron-self-Compton (short-dashed), 
and external Compton (dot-dashed). 

Inverse-Compton peak 

Accounting for attenuation of the gamma rays due to pair pro- 
duction in collisions with photons from the integrated light of 



galaxies, brings the emitted spectrum into agreement with mod- 
els assuming inverse-Compton scattering of external photons as 
the dominant radiation process. The emitted spectrum also ap- 
pears to be very similar to the spectra of many high-peaked BL- 
Lac objects, although their corresponding synchrotron peak is 
at much higher energies. In the framework of inverse-Compton 
scattering models this implies that external photon fields may 
act as scattering target, although seemingly not to the extent to 
create internal opacity due to pair production in photon-photon 
collisions W7\ . The inverse-Compton models predict a steep 
decline for two reasons: for one, the spectrum of the highest 
energy electrons, judging from the synchrotron spectrum, tends 
to be steep; and second, VHE electrons scatter with ambient 
photons in the Klein-Nishina regime, reducing the scattering 
efficiency. An EBL larger than the integrated light of known 
galaxies would imply that the generic leptonic emission mod- 
els must undergo a major revision, perhaps involving proton 
synchrotron and cascade emission. The latter would shed new 
light on the correlation between active galactic nuclei and the 
directions of the highest energy cosmic rays ( |?5] l. 
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